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Abstract 

The adsorption of potassium on the Cu(lll) surface in a (2 x 2) pattern has 
been simulated with all-electron full-potential density functional calculations. 
The top site is found to be the preferred adsorption site, with the other highly 
symmetric adsorption sites being nearly degenerate. The bond length from 
potassium to the nearest copper atom is computed to be 2.83 A. Population 
analysis and density of states indicate that there is no evidence for covalent 
bonding so that the binding mechanism appears to be a metallic bond. 
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I. INTRODUCTION 



The study of adsorbates on metal surfaces is of high importance for scientists as well as 
for industry with potential applications, for example in catalysis. Alkali metals are relatively 
simple systems to consider as adsorbates. Initially, it was assumed that they would always 
occupy high-coordination sites. It was therefore a big surprise when a first system was 
discovered where Cs, adsorbed on the Cu(lll) surface, occupied the top site and not the 
threefold hollow site0. Meanwhile, further systems of alkali metals with top adsorption have 
been found — typically heavier alkali atoms (from K on) on closed packed surfaces (for a 
review, see references The situation is thus more complex than in the case of halogens 
where the adsorption site on closed-packed surfaces is usually the threefold hollow site. 

Due to the increase in computational power, the theoretical study with methods based 
on density functional theory (DFT) of these adsorbate systems has become feasible in the 
last few years. For example, systems such as Na and K on the Al(lll) surface! or Na on 
the Cu(lll) surface! have been studied recently. 

In this article, we study the adsorption of potassium on the Cu(lll) surface as a model 
system for alkali adsorption on metal surfaces. This is an extension of earlier work of 
halogens on metallic surfaces (CI on Cu(lll)i and CI on Ag(lll)i). For the system Cu(lll) 
(2x2)-K, adsorption on the top site was observed experimentallyil. The addressed questions 
are therefore whether first principles simulation is able to reproduce the preferred adsorption 
site, the magnitude of the energy splitting of different highly symmetric adsorption sites, 
the geometry, the charge of the adsorbate, and the mechanism of the bond. 

The applied methods are all-electron full-potential calculations with a local Gaussian 
basis set. Mainly gradient corrected density functional calculations were performed, the 
energetically most favorable structure was additionally optimized with the local density 
approximation, for comparison. 
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II. COMPUTATIONAL PARAMETERS 



A local basis set formalism was used where the basis functions are Gaussian type orbitals 
centered at the atoms as implemented in the code CRYSTAlE. For Cu, a [6s5p2(i] all-electron 
basis set as used in previous studies on copper metal and chlorine adsorption on copper metal 
was employed^!. The K all-electron basis set from work on alkali halides was chosen0, with 
the two outermost sp-exponents replaced by two sp-exponents optimized for the free K atom 
(0.29 and 0.08 a.u.) and an additional (/-exponent (0.50 a.u.), resulting in a [5s4pl(i] basis 
set. Compared to the calculations on alkali halides, more diffuse basis functions are necessary 
because the potassium atom is not fully ionized and a description of its 4s and 4p states is 
required. A further enlargement of the potassium basis set is achieved indirectly because 
basis functions from neighboring centers are used to describe the charge distribution. Ill The 
bulk of the calculations was done with the gradient corrected exchange and correlation 
functional of Per dew and Wang (GGA)§. For the energetically most favorable site, the 
geometry was re-optimized at the level of the local density approximation (LDA) with Dirac- 
Slater exchang and the Perdew-Zunger correlation functionallii, for comparison. 

An auxiliary basis set was used to fit the exchange-correlation potential. For copper, it 
consisted of 12 even tempered (i.e. the ratio from one to the next higher exponent was kept 
fixed) s-exponents and 12 p-exponents in the range from 0.1 to 2000, 8 (i-exponents and 8 
/-exponents in the range from 0.1 to 100. For potassium, the same s- and p-exponents were 
chosen, and 5 (i-exponents in the range from 0.8 to 100. 

The adsorption was modeled by using a three-layer copper slab with the copper atoms 
arranged as in the face-centered cubic (fee) lattice, at the GGA-optimized bulk Cu lattice 
constantil of 3.63 A. Potassium was adsorbed on one side of this slab. This is presently the 
computational limit for desktop workstations with this approach. This is similar to earlier 
studiesi'i where three layers of metal and a single-sided adsorbate layer were found to be 
sufficient for a theoretical description. A supercell approach with a (2x2) adsorption pattern 
as in the experiment was usedi. This slab was not repeated in the third dimension so that 
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the model is truly two-dimensional. The copper atoms in the top layer were allowed to relax 
vertically, an additional horizontal shift as observed in the adsorption of K on Ni(lll)lll 




would, however, require too much numerical effort to be simulated. For comparison, the 
vertical relaxation of the copper atoms was simulated in two different ways: simulations 
were performed where only a uniform relaxation of the top copper layer was possible, and 
simulations where a different vertical relaxation for the copper atoms in the top layer was 
possible (i.e. substrate rumpling). 

Four adsorption sites were considered (see figure |l|): the top adsorption site with K 
sitting vertically above a copper atom in the top layer, the bridge site with K sitting above 
the middle of two copper atoms in the top layer, and two different threefold hollow sites 
where the potassium atoms are placed vertically above a copper atom in the second (third) 
copper layer. These threefold hollow sites can not be distinguished when only the adsorbate 
and the positions of the atoms in the top copper layer are given, but instead the position of 
the atoms in the second and third layer under the adsorbate layer must be considered. The 
first of the threefold hollow sites is referred to as the hep (hexagonal close-packed) hollow 
site because the potassium atom is vertically above a copper atom in the second copper 
layer. The second of the threefold hollow sites is referred to as fee hollow site because the 
potassium atom is vertically above third layer copper atoms. In the following, we will use 
the notion hep (fee) hollow site to distinguish between these two adsorption sites. This is 
closely connected to the notion of hep (versus fee) lattices where there are atoms two (versus 
three) layers vertically below atoms in a certain layer of a solid. 

With a Pack-Monkhorst net of shrinking factor 16, the number of points in the reciprocal 
lattice was 30 for the fee hollow site and the hep hollow site, and 73 for the bridge and top site. 
The Fermi function was smeared with a temperature of 0.01 (li?/i=27.2114 eV) to make 
the integration numerically more stable. For an extensive test of the various computational 
parameters for metallic systems with a local basis set technique, see reference |1^. 
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III. SUMMARY OF PREVIOUS EXPERIMENTAL AND THEORETICAL 

RESULTS 



The adsorption of K on Cu(lll) has been studied with surface-extended X-ray absorption 
fine-structure (SEXAFS) measurementsi. At a nominal coverage of 0.25± 0.01 monolayer, 
a (2x2) pattern was observed. The adsorption site was identified as the top site and a 
bond-length from K to the nearest Cu atom of 3.05±0.02 A was deduced. 

Cluster models^ with Cur as the adsorbate gave in Hartree-Fock and density functional 
calculations the top site as the preferential one. However, the energy splitting was found to 
be relatively large (e.g. ~ 0.5 eV at the DFT level), the computed bond length was in good 
agreement with the experimentally observed. The computed K charge was relatively large 
with 0.7 - 0.8 |e| at the DFT level. Similar results for the energy splitting and bond length 
were then found when potentials were fitted to model the interaction of the adsorbate with 
the infinite surface. 

GGA calculations have been performed for the systems Cu(lll) (2x2)-Na and Cu(lll) 
(3/2x3/2)-Nai. It was found that hollow sites were preferred. Substrate rumpling was found 
to be strongest for the top site; in general copper atoms below sodium atoms moved into 
the substrate, whereas copper atoms without sodium in the vicinity moved outwards. The 
top site was the one which gained most energy from rumpling and thus the energy difference 
between the different adsorption sites was reduced. 

The adsorption of sodium and potassium on the Al(lll) surface has been studied theo- 
retically with the ldaI and it was found that, whereas sodium adsorbed on a substitutional 
site, the fee hollow site and the top adsorption site were degenerate for potassium for both 
structures considered (-\/3 x ^/3i?30° and 2x2). Again, rumpling was most important for 
the top site. 
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IV. RESULTS FOR K ADSORPTION 



In this section, the results of the calculations on the system Cu(lll) (2x2)-K are pre- 
sented and discussed. Firstly, in tables |I| and y, the adsorption geometry and energy are 
displayed. We notice that the adsorption at the top site is energetically most favorable, i.e. 
the experimental situation is reproduced. For this site, the computed adsorption energy per 
potassium atom is 47 niEh (1.28 eV). However, the other considered sites are only slightly 
higher in energy (by less than 1 mEh). This is close to the limit of what can be resolved and 
there are still parameters which might have an influence on the prediction of the adsorption 
site (e.g. numbers of layers in the slab model, a possible horizontal relaxation of the copper 
atoms, basis set deficiencies due to the lack of availability of /-functions). The magnitude 
of the energy splitting is thus much lower as previously found for halogendli, but in good 
agreement with findings for Na on Cu(lll)i and Na or K on Al(lll)i. 

The value for the adsorption energy is in the range of that for potassium on graphite 
(0.48 to 0.89 eV, depending on the coverage and the relaxationEl) , or that for Na on Cu(lll) 
which was computed to be 1.8 to 1.9 evi. 

The binding energy of a free K layer (with a K-K nearest neighbor distance as in the 
adsorbate system, 5.13 A) was computed for comparison and a value of 0.0229 E^, i.e. 
0.62 eV was obtained. This is in reasonable agreement with previous calculations (0.56 
eVi, LDA, K-K nearest neig hbor distance 4.94 A; or 0.83 evS, LDA, K-K nearest neighbor 
distance 4.52 A; in the latter work the energy of the free atom was obtained from a calculation 
without spin-polarization). 

We notice that the K-Cu bond length increases from 2.83 A (top site) < 3.04 A (bridge 
site) < 3.11 A (fee, hep hollow site). This is in agreement with the argument that the bond 
should be the stronger and the bond length therefore the shorter, the lower the number 
of copper neighbors^. It is consistent with results from simulations for halogens on metal 
surfacesi'l. In addition, we note that those atoms in the Cu surface, which are closest to the 
K atom, relax inwards. This relaxation is strongest for the structure with the lowest number 
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of nearest neighbors underneath the K atom — in full agreement with the findings in the 
case of Na on Cu(lll)l. The order of magnitude (0.165 A inwards relaxation of copper atom 
1 under the potassium atom, relative to atoms 2,3,4) is in line with Na on Cu(lll) (0.1 A)eI 
and with one experiment for K on Ni (111) (0.12 A)0, less with the other experiment for K 
on Ni(lll) (0.01 ± 0.09 A)!!. 

The computed bond length is shorter than measured in the experiment (3.05 ± 0.02 
A)i. To estimate the errors due to the choice of the functional, calculations with the LDA 
were performed for the top site, for comparison. The LDA bond length (2.73 A) is slightly 
shorter as the GGA bond length. If we assume a copper radius of 3.63/2v^ A, we obtain 
an effective radius, at the GGA level, of 1.55 A which is in the range of the experimentally 
observed value for top adsorption of potassium (1.77 A for K on Cu(lll), 1.57 ... 1.77 A for 
K on Ni(lll)) — see reference ^ and references therein. 

In addition, calculations for the same adsorption sites were performed when rumpling 
was not allowed, but instead only a uniform relaxation of the copper atoms in the top Cu 
layer (table |T|). In these simulations, the fee hollow site was found to be lowest in energy, 
with the hep hollow site and the bridge site of the order of ~ 0.2 mEh higher in energy, i.e. 
nearly degenerate. The top site was highest in energy, but only by less than ImEh, however. 
The bond lengths are relatively similar to the case where rumpling was allowed, the atoms 
in the top Cu layer relaxed inwards by ~ 0.03 A. We can thus conclude that rumpling is 
crucial to obtain the top site as the preferred site — the top site gains > 2 mEh due to 
rumpling, whereas all the other sites gain less than 1 mEh. 



In table |Tl|, MuUiken charges for the potassium atom are given, projected on the different 
basis functions. The free atom in its "^S state has 7 s and 12 p electrons, obviously p^, Py 
and Pz are degenerate. In the free and neutral K layer, p^, py and Pz population slightly 
increase whereas the s population decreases. Finally, in the adsorbate systems, p^, Py and 
Pz charges are slightly above 4, the number of s-electrons is ~ 6.5, and the total charge is 
~ +0.1 |e|. Charges in p^ and Py orbitals are identical except for the bridge site - as in the 
case of CI on Ag(lll)i, the Px orbital which has more overlap with the substrate layer (in 



our choice of geometry — figure |T]) is slightly less occupied than the Py orbital. However, the 
difference is much smaller than in the case of CI on Ag(lll), and it is hardly visible in the 
density of states (DOS) — see the discussion in the following paragraph. In addition, the 
overlap population is very small: for K and next neighbor Cu, it is 0.06 |e|, and it rapidly 
decreases for the copper atoms further apart. All this indicates that there is no evidence for 
hybridisation and covalent bonding. 

In figure |], the DOS, projected on the potassium basis functions, is displayed for the top 
adsorption site. The Fermi energy is in a regime with a relatively large contribution from K 
4s and 4p bands, so that the overlayer is clearly metallic. It is interesting that the projected 
DOS does not depend on the adsorption site: it looks virtually identical for top, hep hollow, 
fee hollow and bridge site. This is different from the case of chlorine as an adsorbate, where 
the projected DOS clearly depended on the adsorption site!. Another feature which was 
found for the bridge site in the chlorine adsorption was the different projected DOS for 
and Py orbitals: however, this is hardly visible for potassium and the DOS for p^ and Py is 
virtually identical. 

Also, the values of the K 3s level (-1.19 Eh with respect to the Fermi level) and 3p levels 
(-0.59 Eh with respect to the Fermi level) are practically independent of the adsorption 
site. This is again in contrast to chlorine: for chlorine, the core eigenvalue depended on 
the adsorption site, and it was correlated with the charge of the chlorine atom0. For K, 
the variation of the Mulliken charges with adsorption site is much smaller (table p|) as 
for chlorine. This is consistent with the finding that the core eigenvalue does not vary for 
the different sites. The Mulliken charge for K at the top site is ~ +0.12, i.e. K is slightly 
positive charged. This is in reasonable agreement with the findings for K on graphite, where, 
depending on the coverage and the relaxation, charges in the range from 0.17 to 0.38 |e|, 
were obtainedlii. 
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V. SUMMARY 



In this article, all-electron full-potential density functional calculations were applied to 
compute properties of the system Cu(lll) (2x2)- K. The top site was found as the preferable 
adsorption site, with the other considered sites being nearly degenerate. Substrate rumpling 
was crucial to obtain this result. The bond length was computed to be 2.83 A, the binding 
energy 1.28 eV. The potassium charge was computed to be -1-0.1 |e|. Charge and projected 
density of states varied only very little for the different adsorption sites, and similarly the 
position of the K 3s and 3p core eigenvalues did not vary. No evidence for a covalent 
contribution to the binding was found, so that the mechanism should rather be a metallic 
bond. 
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TABLES 

TABLE L Adsorption of K on the Cu(lll) surface, GGA results, if not stated otherwise, di 
is the relaxation of those atoms in the Cu top layer which are nearest neighbor to the adsorbed 
potassium, d2 is the relaxation of those atoms which are not nearest neighbor to the adsorbed 
potassium. The relaxations are given with respect to the bulk value, a negative sign means a 
contraction relative to the bulk. The bond length cZk-Cu nn is the distance from potassium to the 
nearest copper atom. The adsorption energy is the difference at Cu{iii) ~ -E'cu{iii) — E^. The 
numeration of the atoms is like in figure 111 



site nearest relaxation atoms relaxation fix-Cu nn EadsorpU, 

atoms to K di [A] not under K ^2 [A] [A] -r^EWI 



hep hollow" 


2,3,4 


-0.045 


1 


+0.06 


3.11 


0.0464 


fee hollow^ 


2,3,4 


-0.050 


1 


+0.07 


3.11 


0.0465 


bridge 


3,4 


-0.068 


1,2 


+0.03 


3.04 


0.0465 


top 


1 


-0.14 


2,3,4 


+0.025 


2.83 


0.0473 


top (LDA) 


1 


-0.13 


2,3,4 


+0.012 


2.73 


0.0660 



" note: Cu atom in third layer vertically under K atom 
^ note: Cu atom in second layer vertically under K atom 
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TABLE II. Adsorption of K on the Cu(lll) surface, GGA results, without rumphng. The 
relaxation of the atoms in the Cu top layer with respect to the bulk is given, the negative sign 
means a contraction relative to the bulk, i.e. an inwards relaxation. The bond length c?k— Cu nn is 
the distance from potassium to the nearest copper atom. The adsorption energy is the difference 

Ek at Cu(lll) - -Scu(lll) - Ek- 



site 
hep hollow 
fee hollow 
bridge 

top 



relaxation [A] 
-0.026 
-0.028 
-0.026 
-0.028 



<^K-Cu nn [M 

3.10 
3.11 
3.02 
2.80 



E, 



adsorption 



E. 



K atom 

0.0457 
0.0459 
0.0457 
0.0451 



free K layer 



■^cohesive — 0.0229 
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TABLE III. Orbital-projeeted ehar 


•ge of K on different adsorption sites. 




site 






eharge, in |e| 








s 


Px 


Py Pz 


d 


total 








with rumpling 






top 


6.484 


4.155 


4.155 4.072 


0.012 


18.879 


hep hollow 


6.481 


4.149 


4.149 4.072 


0.012 


18.863 


fee hollow 


6.481 


4.150 


4.150 4.072 


0.012 


18.865 


bridge 


6.481 


4.147 


4.155 4.072 


0.013 


18.867 








without rumpling 






top 


6.493 


4.151 


4.152 4.079 


0.013 


18.886 


hep hollow 


6.483 


4.149 


4.149 4.072 


0.012 


18.863 


fee hollow 


6.484 


4.149 


4.149 4.074 


0.012 


18.866 


bridge 


6.485 


4.146 


4.151 4.074 


0.014 


18.868 








free K layer 








6.806 


4.094 


4.094 4.006 


0.000 


19.000 








free K atom 








7.000 


4.000 


4.000 4.000 


0.000 


19.000 
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FIGURES 

FIG. 1. The structures considered for K, adsorbed on the Cu(lll) surface, at a coverage of 
one fourth of a monolayer, (2 x 2) unit cell. The copper atoms in the top layer are displayed by 
open circles. In the top adsorption site, potassium occupies the sites above the copper atoms with 
number 1 (filled circles). Alternatively considered adsorption sites were the threefold hollow sites 
above atoms 2,3,4 (fee or hep hollow, circles with horizontal lines) or the bridge site above atoms 
3 and 4 (circles with horizontal and vertical lines). Note that fee and hep hollow sites can not be 
distinguished in this figure. The positions of the second and third layer copper atoms underneath 
the potassium layer are necessary to differentiate between fee and hep hollow sites. 
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FIG. 2. Density of states, on K projected, top adsorption site. The displayed structure is the 
broad K Asp band. The density of states, projected on all K basis functions, is shown with a full 
line. A projection onto the K s-basis functions only is shown with a dashed-dotted line. The Fermi 
energy is indicated with a dotted line. 
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